Slow shocks and conduction fronts from Petschek reconnection of 
skewed magnetic fields: two-fiuid effects 



D.W. Longcope/'^ and S.J. Bradshaw^'*^ 

1. University of Maryland, College Park, MD 20742 
2. Permanent address: Department of Physics, Montana State University, 

Bozeman, MT 59717 
3. NASA Goddard Space Center Flight, Creenbelt, MD 20771 
4- Department of Physics and Astronomy, Rice University, Houston, TX 77005 

Draft: June 3, 2010 

ABSTRACT 

In models of fast magnetic reconnection, flux transfer occurs within a small 
portion of a current sheet triggering stored magnetic energy to be thermalized 
by shocks. When the initial current sheet separates magnetic fields which are 
not perfectly anti-parallel, i.e. they are skewed, magnetic energy is first converted 
to bulk kinetic energy and then thermalized in slow magnetosonic shocks. We 
show that the latter resemble parallel shocks or hydrodynamic shocks for all skew 
angles except those very near the anti-parallel limit. As for parallel shocks, the 
structures of reconnection-driven slow shocks are best studied using two-fiuid 
equations in which ions and electrons have independent temperature. Time- 
dependent solutions of these equations can be used to predict and understand 
the shocks from reconnection of skewed magnetic fields. The results differ from 
those found using a single-fluid model such as magnetohydrodynamics. In the 
two-fiuid model electrons are heated indirectly and thus carry a heat flux always 
well below the free-streaming limit. The viscous stress of the ions is, however, 
typically near the fiuid-treatable limit. We find that for a wide range of skew 
angles and small plasma (3 an electron conduction front extends ahead of the slow 
shock but remains within the outflow jet. In such cases conduction will play a 
more limited role in driving chromospheric evaporation than has been predicted 
based on single-fluid, anti-parallel models. 



Subject headings: MHD — Sun: flares — Sun: magnetic fields 
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Introduction 



Magnetic reconnection was proposed to explain the conversion of magnetic energy to 
other forms i r i solar flare s (iGiovanellil Il947l ). In an early detailed model of reconnection 
( iParkerl 119571 : ISweetl Il958l ) plasma was accelerated and heated all within the same current 
sheet that initially stored the magnetic energy. Reconnection i n this mode wa s, however, 
deemed too slow to explain flares. In an alternative proposed by iPetschekl ( 1l964j ) the recon- 
nection process decomposed the initial current sheet into a series of waves and shocks. Flux 
transfer occurred within a very small portion of the original sheet, the diffusion region, while 
acceleration and heating occurs in the much larger slow magnetosonic shocks. 

While it was initially analyzed in steady-state, subsequent work has shown the Petschek 
scenario applies more generally. Even a transient transfer of magnetic flux across a pre- 
existing current sheet will, if it occurs within a localized portion, produce the same set of 
shocks proposed by Petschek (ISemenov et al.l Il983l : iBiernat et al.l 119871 : iHeyn fc Semenov 
19961 : lErkaev et al.ll2000l : iNitta et al.ll200ll ). When the flux transfer is impulsive, the current 
sheet decomposes progressively as the shocks propagate away from the flux transfer point. 

In solar flare models, the heat from a slow shock is expected to b e be conducted along 
fleld lines ahead of the shock itself (IForbes et al.lll989l : lTsunetalll996l ). Except for possible 
non-thermal particles ( iKarpen fc Devord 119871 ) . this conduction front (CP) was predicted 
to be the flrst effect from reconnection energy release to reach the l ower atmosphere and 
drive chromospheric evaporation. A conduction front was identifled by lYokoyama fc Shibata 
( 11997) in numerical simulations and later demonstrated to be an effective driver of evapora- 
tion ( lYokoyama &: Shibatalll998(). The structure o f the CP in this particular simulation has 
been recently analyzed by ISeaton fc Porbej (120091 ) using a boundary layer method. 



Conduction fronts precedi ng shocks are well-known features in hydrodynamics (IBeckei 



1922 



ThomaslEgii : lGradlll95ll ) and magnetohydrodynamics (MHD. lKennellll988l : IXu fc Porbed 



19921 ). when thermal conductivity greatly exceeds viscosity. In an unmagnetized plasma vis- 
cous stress arises from ion momentum while heat flux is carried primarily by the faster-moving 
electrons. The same is true of transport parallel to the fleld in a magnetized plasma. In ei- 
ther case, the ratio of viscosity to thermal conductivity, the Prandtl number, is proportional 
to the ratio of the thermal velocities of the two species, which is the square-root of their 
mass ratio and hence very small. 

The reasoning above assumes, as does any single- fluid hydrodynamics such as MHD, 
that the temperatures of both species are equal everywhere and always. More sophisticated 
two-fluid treatments, where the electrons and ions both have nearly Maxwellian velocity 
distributions but with separate temperatures, show the equal-temperature assumption to be 
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inval id throughout a typical parallel shock ( iJukea 119571 : iTidmanI Il958l : iJaffrin fc Probstein 
1964 ). Solutions do include features similar to those of single- fluid models, although modified 
by the disparate temperatures. There is an inner sub-shock, also called an ion shock, where 
ion kinetic energy is thermalized and a preceding CF where heat is carried by electrons. 
In two-fluid models, however, heat must be transferred from ions to electrons at a collision 
frequency far smaller than those responsible for either viscosity or conduction. This transfer 
occurs downstream of the two-fluid shock, in a region called the ion cooling region. Owing 
to the lower frequency of exchange this region is far larger than the other shock regions; it 
is necessarily absent from single-fluid treatments. 

Given the critical role played in flares by shocks and their associated CFs, it is essential 
to model them using a two-fluid formalism. The two-fluid effects will depend on the magnetic 
field at the slow magnetosonic shock outside the flux transfer region. This will depend in 
turn on the field initially separated by the current sheet — the field lines joined together by 
the reconnection. Petschek's original model applied to a current sheet separating perfectly 
anti-parallel fields and had two "switch-off" shocks across which the field direction and 
strength changed dramatically. The model was subsequently exten ded to skewed cases where 
the initial cu r rent sheet separated fields differing by ^ 180° (jPetschek fc Thornd 119671 : 
SowardI Il982l : ISkender et al.ll2003l ). In this case each slow magnetosonic shock is replaced 
by two different shocks: a rotational discontinuity, where the field direction changes at 
constant magnitude and temperature, and a slow shock where kinetic and magnetic energy 
are partially converted to heat. The CF will originate from the slow shock so this is where 
we expect ion and electron temperatures to differ in a two-fluid solution. 

Most studies of reconnection-generated CFs have used single-fluid equations (MHD) 
and worked with the special case of anti-parallel reconnection {A6 = 180°). The present 
study examines the structure of CFs for the more general case, A^^ 7^ 180°, using two-fluid 
equations. It turns out that for most such cases, say A6 < 120°, the slow magnetosonic shock 
is very similar to a parallel shock or an unmagnetized hydrodynamic shock. We therefore 
perform our analysis using two-fluid eq uation s for an unmagnetized p las ma and find steady- 
state solutions very similar to those of I Juked (Il957l ). iTidmanI ( 119581 ) or iJaffrin fc Probstein 
( 119641 ). To apply these results to Petschek reconnection, however, it is necessary to use time- 
dependent solutions. We find that such evolution is dominated by the very slow thermal 
equilibration between electrons and ions. In both respects the shocks differ from single-fluid 
solutions. 



-4- 



2. Petschek reconnection and hydrodynamic shocks 

Reconnection occurs across a current sheet (tangential discontinuity) separating field 
lines of equal magnitude, Bq, but differing in direction by A^. The basic premise of the 
Petschek model is that a localized electric field occurs within an otherwise steady current 
sheet. Original versions of the model attempted to attribute this electric field to Ohmic 
resistivity, but it has subsequently been found that any mechanism producing E with a 
component along the magnetic field will generate very similar external structure — structure 
outside the "diffusion region" designated X in Fig. [H That external structure consists of 
several shocks emanating from the diffusion region. In a strictly steady-state interpretation 
these are standing shocks, however, they are more realistically understood as the ends of 
propagating shocks encompassing the fi eld lines which were topologically changed by the 



inner electric field (ISemenov et al.lll983[ ). The diffusion region can affect very little of the 
plasma directly, owing to its very small size, so it is the shocks launched by the reconnection 
which convert the magnetic energy stored in the current sheet to other forms. 

Two types of MHD shocks are launched from the diffusion region, r otational discontinu- 



ities ( RDs) and slo w magnetosonic shocks (SSs) also called slow shocks (jPetschek fc Thorne 



19671 : ISowardlll982l ) . Magnetic tension accelerates the plasma at the RDs creating an Alfven- 
speed outflow jet. Slow shocks occur within this jet, heating and compressing the plasma. 
Thus only the inner portion of the outflow jet is hot. In the special case of anti-parallel 
flelds (A^^ = 180°) these two distinct shocks merge into a single slow shock of the switch-off 
type, often associated with the Petschek model. Naturally one cannot assign, in this singular 
limit, different roles to the different shocks, so one often hears that slow shocks accelerates 
the outflow jet as well as heat the plasma. In the general case, however, slow shocks only 
heat while RDs accelerate the jet. 

There is a relatively small inward velocity, the "inflow" which we call vq, in the outer 
region (regions and 0' of Fig. [1]). In genuine steady state models this is often interpreted 
as an external driver of reconnection, however, in transient or spontaneous model s it is part 



of a fast magnetosonic rarefaction wave launched at the same time as the shocks (ILin fc Lee 



1994J : iNitta et al.ll200ll ). In that case the flow is caused by the reconnection and does not 
cause the reconnection. While it is a relatively slow flow [vq ^ ^a,o) it has been the focus of 
intensive study since reconnection models were flrst proposed — indeed, how it compares to 
VAfi is generally regarded as the crux of the fast reconnection problem. In spite of this focus, 
many aspects of the post reconnection dynamics, including heating and energy release, are 
largely independent of vo- We focus exclusively on these energetics, occurring in regions 1 
and 2, and hereafter ignore any effects strictly dependent on the inflow vq. 



The structures of the different regions follow from the standard jump conditions of 
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Fig. 1. — Petschek reconnection of skewed magnetic fields. Reconnection occurs within a 
portion of a large-scale current sheet (far left). Two perspectives of the reconnection region 
are shown: the end-on view (center, whose horizontal scale is exaggerated for clarity) and 
the face-on view (right). Flux transfer occurs along a line X (grey) along the ignorable 
direction (z). Shocks, RD', SS', SS and RD (thick dashed lines) originate there, dividing the 
reconnection region into sub- regions labeled 0', 1', 2, 1 and (inside boxes). The central 
region, 2 (shaded grey) contains plasma which has been heated and compressed by the slow 
shocks. Reconnected field lines (thin sohd lines) bend at each shock. A single fiux tube is 
called out in grey in each panel. It passes though the shocks at points R', S', S and R, left 
to right. 
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the MHD shock types. The RD is a non-hnear manifestation of a shear Alfven wave and 
thermodynamic quant ities, density, pressure and temperature are all continuous across it 



( iPriest fc Forbes! l2000l ). Thus the RDs convert magnetic energy to bulk kinetic energy. The 
magnetic field changes direction but not magnitude there so the magnetic energy release is 
due to field line shortening alone. The slow shocks are the site of all heating; it is there 
that the temperature and density both jump to higher post-shock values, T2 and p2- The 
magnitude of the field weakens across slow shock so B2 < Bi. Figure [2] shows that for typical 
coronal conditions, /3o = 0.01, the amount of weakening is modest except for the anti-parallel 
case where the shock can be considered a "switch-off" shock. The RD rotates the magnetic 
field to enhance the ignorable component: B^i > B^q. The slow shock then weakens this 
enhanced component, but generally leaves it greater than it started: B^i > Bz2 > Bzq. 

With only a modest weakening of the magnetic field, the energy for plasma heating 
comes largely from bulk flow within region 1. The RDs create a flow component v\\^i along 
the local magnetic field. To see this consider the RD from a frame moving with point R in 
Fig. [U This point moves along the unperturbed magnetic field at speed f a,o- In this co- 
moving frame the plasma flows parallel to the field line at constant speed changing direction 
by the angle f2o,i — the angle between Bq and Bi. Neglecting the small external inflow 
(vq — 0) the flow speed in the co-moving frame is that of the reference frame, fA,0; on both 
sides of the RD. Returning to the lab frame we find the parallel component in region 1 to be 



^11,1 = 



^^A,o - wa,o cos (^0,1 = 2t;A,osin2(fio,i/2) 



While the Lorentz force cannot accelerate parallel to the local field, the net effect of a bend 
is acceleration along its bisector. This produces flow with a component parallel to the field 
on the non- stationary side. Indeed, such inward parallel flow is necessary for shortening the 
magnetic field lines. 

When the external plasma is characterized by very small Pq, the parallel flow speed is 
typically greater than the sound speed. 



M, = ^ = 2^sin2((]o,i/2) = J^sin2((]o,i/2) , (2) 

Cs,l Csfi V 7Po 

where 7 = 5/3 is the ratio of specific heats. When /3o ^ 1 each RD produces a flow at 
or above Mach 1, directed toward the center. More than any magnetic change, it is the 
colhsion of these supersonic flows that accounts for the slow shocks. Indeed, Fig. |2] shows 
the magnetic field changes direction very little across these shocks so they are akin to parallel 
shocks. Except for nearly anti-parallel cases AO ~ 180°, the SSs are essentially hydrodynamic 
shocks from pistons moving Mj times the ambient sound speed. The top axis in Fig. |2] shows 
Mj corresponding to the current sheet angle, A6, and the dashed curve is the temperature 
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Fig. 2. — Jumps across the slow shock as a function of reconnection angle for (3o = 0.01. 
The bottom panel shows the jump of the magnetic field magnitude, B2/B1 on the l eft axis, 
and d irection, f2i 2 on the right axis. The s olid curve is from the 2.5d stea dy model of lSoward 



(Il982l ). as reported by lForbes et al.l ( Il989l ) and IVrsna. 



are from a Id Riemann problem solution of iLin fc Lee 



c &: Skenderl (120051 ) . and the diamonds 



(119941 ). The dashed curve is the angle. 



read off the right axis. The top panel shows the temperature ratio for the same two models 
(solid and diamonds). The parallel flow speed, f|| 1 for the Id Riemann solution is used to 
compute a hydrodynamic mach number, Mj, plotted along the top axis. The corresponding 
temperature ratio from a Id hydrodynamics shock, is plotted as a dashed curve. 



ratio for a one-dimensional hydrodynamic shock at that Mach number (jCourant Sz Friedrichs 
1948h . 



The jump conditions across MHD shocks are dictated by conservation of mass, mo- 
mentum and energy through the Rankine-Hugoniot conditions. Transport effects, such as 
viscosity and conductivity, conserve all of these and therefore do not modify jump conditions. 
They do, however, determine the "inner struct ure" of the jumps including the length scale 
over which the quantities change (jKennellll988l ). The largest of the transport coefficients is 
thermal conduction parallel to the magnetic field. This conducts heat from region 2 into 
pre-shock regions 1' and 1, creating a conduction front (CF) ahead of the slow shocks as 
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shown in Fig. |3l It is not immediately clear under what conditions these fronts might extend 
far enough to reach beyond the RDs. They clearly must for anti-parallel fields since the 
RD and SS coincide — the only circumstance in which Petschek-driven CFs were previously 
studied. 



face-on view 

end-on view 




Fig. 3. — Expected location of the conduction fronts (CF) ahead of slow shocks. Reconnection 
configuration is the same as Fig. [TJ Additional grey regions represent the CFs and arrows 
within them show the direction of heat flux. An inset along the bottom of the face-on view 
shows the distribution of temperature along the flux tube. It has constant value, T2 within 
region 2 and decrases smoothly to Ti through the CF. The temperature gradient drives the 
heat flux, — kVT, as indicated by the arrows. 

Instead of tackling this problem using the full equations of MHD we make use of the close 
relation, for most skewed fields, between the slow shocks and hydrodynamic shocks. Since 
transport effects do not modify the jump conditions we know that cases where < 120° 
will involve collision along a flux tube bending by not more than 3° across the SS. This is 
so close to being a straight tube that we apply the analysis of a piston in a straight shock 
tube. 

The initial energy in such a system is primarily the bulk kinetic energy of the ions 
composing the super-sonic flow. Its deceleration at the shock converts most of this kinetic 
energy to thermal energy — thermal energy of ions. The ions must then transfer some of 
their heat to the electrons whose much higher conductive flux creates the CFs. It is thus 
apparent that the internal structure of a hydrodynamic shock must be analyzed in terms 
of separate ion and electron fluids. Steady-state analyses of this kind were performed by 
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Jukea (119571 ). iTidmaru (Il958[ ) and iJaffrin fc ProbsteinI ( ll964l ). We repeat this in §4, below, 
adapting it slightly, before analyzing the transient behavior in §5. Owing to the slowness 
with which ion heat can be transferred to electrons, the steady state would be achieved only 
gradually if ever. In §6 we place the one-dimensional, time dependent shock solutions into a 
steady-state Petschek reconnection jet to ascertain the effects of two-fluid equations on slow 
shocks and CPs. 



3. The Shock tube problem 

To model the structure and development of the slow shock we solve for hydrodynamic 
flow into a straight shock tube with a fixed end at s = 0, representing the center of region 
2. The fluid enters the tube from the right (s > 0) with properties of the ambient, pre- 
reconnection plasma, pi and Ti, (none of them were changed at the RD). Although we treat 
the electron and ion fluids separately, we assume they are in equilibrium within the inflow, 
Tel = Til = Ti. The flow enters the tube with a leftward velocity v = —v\\^i flxed by the RD 
according to Eq. ([1]). The flow actually originates from the RD, point R, receding at speed 
va,o- We assume that this outruns any effect from region 2 and take the uniform inflow 
region to extend to s — t- +oo. 

In the view of this shock tube picture we use the term inflow for the inward flow, 
V = —"^11,1, originating from the open end of the shock tube. It should be emphasized that it 
represents a component of vi in Pig. [H which is itself inside what is often termed the outflow 
jet. More signiflcantly, it is not the same as the reconnection inflow, vq, which goes into 
the outflow jet from outside. Nor is our shock-tube Mach number, Mj, related to the Alfven 
Mach number often used in reconnection theory. Hereafter we consider only one-dimensional 
hydrodynamics in a straight tube. 



3.1. Two- fluid equations 

Por the reasons elaborated above we use separate equations for ion and electron fluids. 
We assume the shock thickness will be much larger than the Debeye length so quasi-neutrality 
is appropriate: Ue = Ui = n. The central region is current-free so we assume no current 
crosses the slow shock which means that the charged fluids move together: Vg = Vi = v along 
the tube (i.e. normal to the shock). The common density and velocity satisfy continuity and 
the ion momentum equation 

dn d , . 
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{rriinv^) - —{Pi + Pe) + ^ {-fXi^ ] , (4) 



ot OS OS OS \3 as 

where rrii is the mass of an average ion (assumed to be singly charged) and /ij is the dynamic 
viscosity of the ion fluid. The electron pressure appears in place of the electric field after 
neglecting inertia in the electron momentum equation 

eE. = . (5) 

n OS 

Since the electron momentum equation lacks an electron viscosity contribution that effect 
is absent from the ion momentum equation, Eq. (jl]). Were it in cluded it would have been 



smaller than the ion viscosity by me/rrii <^ 1. It has been found by lJaffrin fc Probsteid (119641 ) 
to have an effect only within a Debeye-length layer of charge separation at the leading edge 
of the CF. 

We assume sufficient collisionality within each species to isotropize their distributions 
leading to scalar pressures, pe = nk-QTe and pi = nk^Ti {k^ is Boltzmann's constant). The 
collisional energy exchange between ions an electrons, at frequency Uei, is slow enough that 
we consider their temperatures to be independent and use two separate energy equations. 
The electron fluid's contribution to bulk kinetic energy is neglected, so its entire volumetric 
energy density is the internal energy |pe- This evolves according to 



^ilP<^) = -§;il^Pe) + v^ + -^^(^K,^ j +lkBnu,,{T,-T,) , (6) 

where is the electron thermal conductivity. The second term on the right hand side (rhs) 
represents work done by the electric field, which has been replaced using ([5]). The third term 
is heat transport by thermal conduction and the final term is heat supplied by the ions. 

The scalar ion pressure evolves according to a similar equation. To facilitate the deriva- 
tion of conservation laws we write the equation governing the evolution of the total ion 
energy, thermal plus kinetic, 

d ,n 1 2\ 9 1 s\ ^Pe d f 4 dv\ 



dt 



d f dTi 

+ Tr i^i^ ) + I^B nUei{Te - Ti) . (7) 



ds V ds 



The second term on the rhs is work from the electric field, with sign opposite the term in Eq. 

since the ions are positively charged. The third term is the work done by viscous stress 
both on the bulk kinetic energy and on the internal energy as irreversible heating. The final 
term is the heating from electrons; this is naturally equal and opposite to the corresponding 
term in Eq. (|6]) since the energy is merely exchanged between the two species. 
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To bring out the fundamental scalings in the shock solutions we rescale all quantities 
using properties of the inflow. The velocity is scaled to sound speed Cgi which we define 
using the single-fluid form 

^2 _ 5 Pel + Pil _ 10 fceTi 
3 ninii 3 

The asymptotic inflow velocity, v —>■ —Mj, as s —>■ +00, is the one free parameter in the 
problem. We rescale temperatures to rriiC^i/kB so the normalized inflow temperature is 
Ti — 0.3. We rescale number density to rii. 

The fundamental length scale in a shock is set by the viscosity. This varies with tem- 
perature so we use the ambient value to define a reference length scale 

3 miTiiCsi VIO Ky V ^ei / 

related to the pre-shock ion mean free path. Times are rescaled to Th — in/cgi which is 
proportional to the ion-ion collision time. 

The local viscosity is proportional to the five-halves power of the local ion temperature, 
which we denote by Xi for notational convenience 

/^i = /^af^j = l^iiXi(Ti) . (10) 

Thermal conductivities have identical temperature dependence which we express using the 
same notation 

J, \ 5/2 



The ratios of viscosity to the different conductivities defines Prandtl numbers 

3 kjiKg^i 

which are independent of temperature. For Spitzer values of these transport coefficients, 
and taking to be the proton mass, these dimensionless quantities arc Pr^ = 0.013 and 
Prj = 0.33. The electron-ion collision frequency has a related temperature dependence, 
which we write as 

4v/2^eVfeT.^-/^,^^^2Ax^_ (13) 



n Smerrii \ nie J 3 riiTi 
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where In A is the Coulomb logarithm. For Spitzer collision rates the coefficient is /i = 0.019. 
The rescaling yields a set of equations including an unchanged continuity equation, (jH]), 



and 



d{nv) 
dt 



d_ 

ds 
d 



d{nT, 



ds 



+ 



d 



d 



+hn\-''\T,~T, 
d_ 
ds 



ds 
d{nT, 



Xe 

Pr. 



dv 
as 

ds 



d 



d 

+ 7^ 



ds \ Pr,- ds 



ds ds 
f,n\-^/\T,-T, 



dv 



(14) 



(15) 



(16) 



These are to be solved with inflow boundary conditions, v = —Mi, Te = Ti = 0.3 and n = 1 
at the right, s — )■ +oo and mirror conditions at s = 0: v = 0, dT^/ds = dTi/ds = dn/ds = 0. 
The only free parameter, Mj, in the inflow boundary condition, represents through Eq. ([2]), 
the Alfven speed and current sheet angle of the reconnection. 



3.2. The single- fluid limit 

It is the small value of the dimensionless constant /i that dictates the slow coupling 
between the electron and ion temperatures. In spite of this fact, it is common to treat shocks 
using single-fluid equations whereby the two temperatures are equated, Tg = Tj = T. In this 
case the momentum equation, f|T^ becomes 



dt ds ds ds \ ds 

where x{T) = (T/0.3)^/^. The net energy equation comes from summing (ITSl) and (fT6ll . 



where the combined Prandlt number, Pr = PrgPrj /(Prj + Prg) = 0.012, includes both elec- 
tron and ion conductivities. The pressure, p = 2nT, is a combination of ion and electron 
contributions now required to be equal. We use these single-fluid equations for comparison 
to a larger body of work on purely hydrodynamic shocks and MHD parallel shocks. 
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4. Stationary shock solutions 

We begin the analysis of Eqs. ([3]) and f|T4|) -f lT6|) by seeking solutions which are stationary 
in a reference frame moving to the right at some speed Vs- Since all the equations are invariant 
under Galilean transformation, the only modifications are to the boundary conditions. The 
upstream conditions, s — )■ +cxd, are the same, = Ti = 0.3 and n = 1, but with a boosted 
velocity 

V ^y-^ = -(^y^+Mi) = -Ml . (19) 

At the left boundary, now at s — t- — oo, we assume all quantities approach uniform values, 
Tg — )■ Te2, Tj — )■ Tj2, n — )■ 77-2 and f — )■ f2. (We do not assume a priori that the post- 
shock electron and ion temperatures are equal but show below that they are.) These post- 
shock values must be found from the solution of the equations. Once found, the velocity 
of the reference frame is then determined by making the downstream velocity vanish in the 
laboratory frame: Vg = |f2| and Mi = Mj + \v2\. 

The solution we seek is time-independent in this reference frame, so = in all 

equations. The continuity equation, Eq. ([3]), is then quickly integrated to yield nv = riiVi = 
—Ml. This can be used to formally eliminate n in favor of v. The momentum equation, Eq. 
(JH]), can also be integrated to yield 

.,Z^,^* , ,20) 

V Ml OS Ml 

after dividing by riiVi and noting that dv/ds in the upstream region. A third conser- 
vation law, conservation of total energy, results from the addition of Eqs. f|T5|) and f|T6|) 



IK + T.) + + ^1 + -^f + . 5r, + iMf (21) 

Equations f l20|) and fl2T]) can be cast as first order equations for the profiles v{s) and 
TAs). The first is 



dv _ Ml 

ds XiV 



T, + Ti + v'+ {Mi + ^^]v 



M 

^[T, + T,-2giv)] , (22) 



after using Ti = 0.3. To cast Eq. ( 1211) as an explicit first order equation for Tj we eliminate 
dv/ds using Eq. f l22p and introduce the heat fiux variable 
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yielding 

^ = -^[l{T^ + T,) + y + 2g{v)-\Ml-l + FjM,] 

= -^^[T, + T,-2%) + |Fe/Mi] . (24) 

For notational convenience we have introduced two quadratic functions of velocity 
g[v) = - {v + Mi + 0.6/Mi) , h{v) = l[v^ + 2(Mi + 0.6/Mi)v + Ml + 3] . (25) 

A final relation must come from either the electron or ion energy alone. Neither of these 
is independently conserved, so the second-order equation can not be simply integrated. The 
electron equation, ( ITSjl . is, however, of first order in heat flux F^. Using Eq. (122|) to once 
again eliminate dv/ds gives the explicit equation 

Equations f l22|) - fl26|) are four first-order equations for the four profiles f (s), Tj(s), Te(s) 
and Ff.{s). Their solution describes the smooth transition from an upstream state, charac- 
terized by = —Ml, Te = Ti = 0.3, and F^ = 0, to the post-shock state down stream. These 
two states must be fixed-point solutions of the equations: solutions for which d/ds = 0. 



4.1. Fixed points 

Fixed point solutions to Eqs. fl22|) - fl26|) are found by requiring d/ds = in each equation. 
According to Eq. fl23|) any fixed point state must have vanishing electron heat flux, F^ = 0. 
Equation fl22l) requires that Tg + Tj = 2g{v) at the fixed point. Using these two in Eq. 
f l26p shows that Tg = Tj at a fixed point solution. Using all of these in Eq. f l2^ shows that 
Te = Tj = h{v). Since Tj = g{v) = h{v), the velocity must be a root of the quadratic 
equation 

h{v) - g{v) = = |(t; + Mi)(4t; + Mi + 3/Mi) . (27) 

There are therefore exactly two fixed point solutions to the coupled set of four differential 
equations (122|) - (!26|) . These are naturally the upstream and downstream state between which 
the shock makes a continuous transition. 

One solution of the quadratic equation, v = —Mi, is the trivial one: the upstream state 
used to set the constants in the first place. The other fixed point solution, 

V2 = -\M, - fMf^ , (28) 
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is the downstream state. The downstream temperature is T2 = g{v2)- The ratios of fixed 
point quantities can be seen to match the well-known Rankine-Hugiot relations for single- 
fiuid hydrodynamics 



V2_ 

Vl 



ni 
n2 
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5Mi 



(29) 
(30) 



Thus it seems that a two-fiuid shock satisfies the same Rankine-Hugoniot conditions as a 
simple hydro dynamic shock. 



From these we find an explicit expression for the infiow Mach number 



Ml 



F2I 



KMi-Mfi) 



(31) 



The dashed line in the top panel of Fig. [2] shows T2/T1, from Eq. ([30 
Mach number, Mj (the top axis). 



versus this infiow 



The internal shock structure is a solution of Eqs. fl22l) - fl26l) diverging away from one 
fixed point and then converging to the other. Finding this solution in practice requires a 
knowledge of the asymptotics of each fixed point solution. The asymptotic behavior can 
be obtained from a linearization of the governing equations about the fixed point. Electron 
and ion temperatures match at each fixed point so the two temperature dependent-functions 
match: Xi = Xe = X- Using this fact we find the linearized equations 
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where all quantities in the matrix are evaluated at the fixed point. The asymptotics near 
fixed point j (j = 1, 2) depends on the four eigenvalues of the matrix which we designated 
— ^4"''* is ascending order. We find that all four eigenvalues are real for both fixed points 
for all Mach numbers we have examined. 



4.2. Stationary solutions 

We find that the upstream fixed point has three negative eigenvalues and a single positive 
eigenvalue, ^4^^ The solution can therefore converge to the upstream fixed point along any 
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direction which combines the three eigenvectors with negative eigenvalues; this direction 
must be orthogonal to the fourth direction, otherwise it would diverge from the fixed point 
as s — 7- +00. Were we to integrate Eqs. (122 p in backwards in s, beginning at the 
upstream point, we would be able to begin along any of the convergent direction; these span 
a three-dimensional sub-space. 

The situation is slightly better constrained at the downstream fixed point at which we 
find two positive and two negative eigenvalues for all Mi > 1. The g-integ ration can be mos t 



easily performed from this point in the manner of a shooting method ( jPress et al.lll992l ). 
Begin at a state displaced from fixed point 2 in a direction mixing, with some angle (j), the 
two eigenvectors with positive eigenvalues. The non-linear system, Eqs. fl22l) - fl26|) . are then 
integrated up to a state in the neighborhood of fixed point 1. The small displacement from 
that fixed point is then contracted with the (left) eigenvector of the diverging direction e^^"* 
to produce a scalar function Z{(f)). In order for the solution to ultimately converge to the 
upstream state we seek a solution orthogonal to 64^^ meaning we want Z{(f)) = 0; this will 
not be true in general. The procedure is iterated, repeatedly integrating from fixed point 2 
with different mixing angles 0, until one is found that satisfies Z{(f)) ~ 0. 

Through this non-linear shooting process it is possib le to find a com plet e solution fo r 



any given Mach number Mi . Following pioneering works of lThomasI ( ll944j ) and iJukej ( 119571 ) , 



we find it useful to integrate in velocity rather than position, writing 

dx 1 dTf, dTg/dx 

dv dv/dx ' dv dv/dx ' 

and so forth. The system is then integrated from v = V2 to v = —Mi. 

Figure H] shows the numerical solution for Mi = 3.5 corresponding to an inflow Mach 
number, Mj = 2.4; it is typical of solutions we find at other Mach numbers. It approaches the 
downstream fixed point (2) along a direction where electron and ion temperature discrepan- 
cies, STe and STi, have opposite signs and are much greater than the velocity perturbation, 
6v. We find this to be very close to the eigenvector with the smaller of the two positive 
eigenvalues, This portion of the shock, where Ti > T^, is called the ion cooling region. 

The two temperatures cross at some point in phase space and then approach the up- 
stream fixed point (1) with Tg > Tj. It is evident from Eq. (I22l) that a monotonic velocity 
profile requires the mean temperature, T = (Te -|- Ti)/2, to be less than g{v). This means 
that the orange curve in Fig. H] must remain below the upper green curve. We find that the 
mean temperature, T approaches g{v) from beneath with the electron temperature relatively 
constant. It then follows T c:^ g{v) as the electron temperature decreases to Ti. This portion 
of the solution is the conduction front, designated CF in the figure. 



-17- 




0.2 0.4 0.6 0.8 1.0 1.2 -800 -600 -400 -200 200 
v/v, s 



Fig. 4. — Stationary shock profile for Mi = 3.5. Tlie bottom panels sliow temperatures for 
electrons (blue), ions (red), their mean (orange) and the single-fluid solution (dashed violet). 
The left panel shows the phase space, T vs. v, both normalized to upstream values. The 
fixed points are shown with asterisks and labeled 1 and 2. Green curves show the quadratics, 
g{v) and h{v) from Eq. fl25l) . The right panel shows the profile against distance, s, in which 
flow is from the right. The top panel shows the velocity (black against right axis — negative 
is up) and density (green against left axis). The profiles are broken into regions which are, 
from right to left, the conduction front (CF), the sub-shock and the ion cooling region (IC). 
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The structure of the solution can be understood in terms of dynamics of the electron and 
ion fluids. Approaching super-sonically from the right the ions slow abruptly to sub-sonic 
in the sub-shock, also called the ion shock. Ion viscosity smooths the jump, so its size is 
governed by fii{T). The kinetic energy of the ions is converted to ion thermal energy through 
viscous heating. Since ion conduction is quite similar to the viscosity (Prj = 0.33) heat flux 
does not broaden the sub-shock appreciably. 

Quasi-neutrality demands that the electrons slow down with the ions, within the sub- 
shock. Due to their much smaller mass, they carry far less kinetic energy than the ions, and 
the sub-shock does not heat them very much. Since only the ions are heated at the sub- 
shock the electrons are cooler than the ions downstream of it. The ions then transfer their 
heat to the electrons over the ion cooling region via the small collision frequency z/gj. Due 
to their very large conductivity the electrons conduct this heat back upstream producing a 
CF of pre-heated electrons in advance of the sub-shock. Its extent reflects the large electron 
conductivity, PrJ^ = 77. The electron pressure in this front partially decelerates the ions in 
advance of the sub-shock, through a polarization electric field. 

Since the solution near fixed point 2 follows the slower of the two positive eigen- 

(2) 

directions, Q , the size and structure of the ion cooling region is dictated by it. In the 
hypothetical absence of all electron-ion collisions, Uei = /i = 0, the downstream state could 
be a fixed point with different ion and electron temperatures provided their mean was g{v2)- 
This indeterminacy gives rise to a zero eigenvalue in the matrix for that hypothetical limit 
(/i = 0). It is evident from inpsection that setting /i = in the matrix renders it singular 
with null eigenvector [0, 1, —1, 0]. Adding this to a solution would change the ion and elec- 
tron temperatures without changing their mean. The null vector is a manifestation of the 
temperature independence in the absence of collisions. 

If the collision rate is now increased from zero, the eigenvalue increases approximately 
with the corresponding term in the matrix 



one another as they converge to their downstream value. Its inverse is the size of the ion 
cooling region Lie ~ ^i/Zi; using the hypersonic (Mi ^ 1) limit of Eqs. fl29|) and f pOj) . 

The other regions in the profile can be analyzed only using the full non-linear solution. 
We define the size of the sub-shock to be twice the separation between the velocity mid-point, 
f = (f 1 + ^2)/2, and the velocity 95% to the downstream value. The extent of the cooling 

(2) 

region is defined to be 3/^3 . All lengths vary significantly with Mach number as shown in 




(33) 



This is the spatial rate at which the electron and ion temperatures exponentially approach 
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Fig. |5l It is evident that the largest of the three regions is the ion coohng region. This is a 
consequence of z/gj being the smallest collision frequency. The sub-shock is consistently the 
smallest of the regions, but is greater than the ambient length, in (i.e. unity), due to both 
the increased viscosity in the hotter plasma and the effects of ion conductivity. 




Fig. 5. — The scales of the three regions in the stationary shock, versus the shock Mach 
number Mi (bottom axis) and inflow Mach number Mj (top axis). The right axis are the 
lengths in Mm when rii = 10^° cm~^ and Ti = 10^ K. Solid curves are, from top to bottom, 
the lengths of the ion cooling region (IC), the conduction front (CF) and the sub-shock. 
Triangles are the values from Fig. H] at Mi = 3.5. Dashed curve is the approximation 
Lie = 0.05M,V/i. 



4.3. Fluid approximation of electron heat flux 



The magnitude of the electron heat flux — KgVTg is given by the dimensionless variable 
Fe times rriinicl i- This classical, local form of heat flux is derived assuming small departures 
from a Maxwellian electron distribution. This app roximation, is ex pected to break down as 
the heat flux approaches a "free-streaming" value (jCampbelllll984j ) 



th,e 



§7 ^nii/me 



niTjTifl^] 



(34) 
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where the quantities inside curly braces are the dimensionless versions. Since both quantities 
scale identically with density and ambient sound speed (the factor in square brackets), their 
ratio will depend only on the rescaled solutions above. Figure M shows that the actual heat 
flux, Fg, is well below the free-streaming limit throughout the shock; in all cases shown 
the peak heat flux is below the limit by a factor less than 2.5%. We thus conclude that a 
classical, local heat flux expression is justified in stationary slow mode shocks and their CFs. 




Fig. 6. — The heat flux in stationary shocks of different Mach number. The top panel shows 
temperature profiles of the electrons (solid) and ions (broken line) for Mi = 3.5 (left), 5.0 
(middle) and 7.0 (right). The horizontal scale is the same for all three. Below these are 
the electron heat flux Fe (solid) and the free-streaming limit, Fe^^ (dashed) for each case. 
Triangles mark the location of peak heat flux. 

The fact that F^ <^ Fe^^-' is equivalent to the statement that the temperature gradient 
scale, £t = \dhiTe/ds\~^ , is much greater than the local electron mean free path, ie, through- 
out the shock. The ratio of the polarization electric field to the Dreicer electric field, the field 
at which thermal electrons would "run away", involves a similar length scale comparison, 
but with the pressure gradient scale in place of It- This is satisfied to the same degree as 
the heat flux, so the polarization electric field is everywhere much less than the Dreicer field. 
These statements are not entirely trivial since the gradient scales in the solution are set by 
diffusive processes and can therefore have length scales similar to the mean free paths which 
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define those processes. 



4.4. The Single-fluid hmit 



Stationary versions of the single-fluid equations can be obtained by setting Tg = Tj = T 
in Eqs. ( 12^ and (^^, and using Eq. (12^ to replace = xidT/ds)/PTe- The coupled 
equations, after some manipulation, are 
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[T - h{v)] 



(35) 
(36) 



with the single-fluid Prandtl number given by the harmonic sum of the electron and ion 
version 1/Pr = l/Fi^ + l/Pr^ = 1/0.012. 

It is evident that these equations have fixed points where T = g{v) = h{v), the same as 
the two-fluid case. The nonlinear system is easier to solve since the upstream fixed point has 
two negative eigenvalues and the downstream fixed point has one positive and one negative 
eigenvalue. One integrates away from the downstream fixed point along the direction of 
the positive's eigenvector (actually anything nearby will converge to that direction). The 
solution then converges automatically to the upstrea m fixed point. Such an a nalysis has 
been presented several times before in neutral ga ses (ThomasI 1944 : Grad 1951), governed 
by ha rd-sphere collisions where xiT) = a/T/Ti. Kennel ( 19881 ) and Guidoni fc Longcope 
( I2OIOI ) used xiT) = (T/Ti)^/^ appropriate to a high-temprature plasma. Such a solution is 
shown along with the two-fluid solutions in Fig. HI 

The single-fluid shock naturally lacks an ion cooling region since the ions and electrons 
have the same temperature by fiat. There is a sub-shock, where the flow decelerates, and 
a CF ahead of it. The phase space curve is roughly horizontal (isothermal) between fixed 
point 2 and the limiting curve T = g{v); this is the sub-shock. The phase-space curve then 
approaches fixed point 1 following T ^ g{v) — the conduction front. 

This simplified two-part structure cannot occur when V2 is to the right of the peak in 
g{v). Setting these two equal gives an equation 



V2 - ^^peak 







(37) 



whose solution is the critical Mach number M^^ i = 3/\/5 = 1.34, below which no CF occurs. 
The same curve, g{v) governs the two-fluid case, so we expect a critical transition there 
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too. The critical shock Mach number corresponds to a very small inflow Mach number 
Mcr,i = 0.44, expected only after reconnection at very small angles: for Pq = 0.01, shown in 
Fig. [2], sub-critical reconnection requires AO < 30°. 

5. Development of the shock 

We expect the post-reconnection flow to approach the stationary solutions found above. 
The inflow boundary of our shock tube is intended to represent the retraction flow established 
in region 1 by the rotational discontinuity (RD). By moving the inflow boundary to s — t- +oo 
we have assumed that the RD recedes rapidly enough that it does not affect the central flow. 
It will not, however, initialize a flow resembling the stationary shock. We propose here a 
simple form for the initialization of our shock tube problem intended to mimic the effect of 
the reconnection which created the tube and the RDs which created the flow. 

We neglect any heating that might have occurred during the reconnection event that 
created the post-reconnection flux tube. We thereby take the thermodynamic quantities to 
be initially uniform: Tg = Tj = Ti = 0.3 and n = ni = 1 in our rescaled variables. This 
approach differs significantly from some previous investigations that focussed exclusively on 
the energetics of the dissipation region — for example on Joule heating. Our approach is, 
however, entirely consistent with the Petschek scenario in which most of the heating oc- 
curs outside the diffusion region in shocks. Moreover, there is sufficient uncertainty about 
reconnection microphysics that no consensus exists on how any initial heat might be appor- 
tioned between electrons and ions. Beginning with uniform ambient properties, as we do, 
is a conservative means of side-stepping the issue altogether. In any event, the diffusion re- 
gion is so small that any initial energy deposited will be quickly overwhelmed by the energy 
thermalized in shocks driven by subsequent retraction. 

For numerical convenience we choose a smooth profile with which to initialize the velocity 

v{s,0) = -Mitanh(s/A) , (38) 

where the initial gradient scale. A, represents the size of the original diffusion region. We find 
that any length smaller than in (i.e. A < 1) leads to similar behavior as the viscosity quickly 
smooths the velocity out to this scale. In the opposite limit, the length scale A defines a 
ballistic time scale, Xb = {X/MijTn — the time on which flow would collapse the profile to a 
singularity. 
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5.1. Time-dependent solutions 

We numerically solve the dynamical equations, ([3]), and f|T4|) -ffT6l) subject to inflow 
boundary conditions at a right boundary, s = L, and the initial conditions described above. 
We take the left boundary to be s = with mirror conditions v{s) = and dn/ds = 
dTf./ds = dTi/ds = there. Alternatively we could have placed a second boundary at 
s = —L from which to drive rightward inflow, v{—L) = +Mi. The center would have been 
mirror symmetric owing to the symmetry of the inflows. 

Figure [7] shows the development of a shock from an inflow of Mj = 2.4 beginning with a 
gradient scale of A = bin. The evolution exhibited here is typical of other Mach numbers. A 
region of high ion temperature develops initially and the region expands into the inflow. The 
central electron temperature begins to rise slowly and gains speed once the ion temperature 
begins falling. It is at this phase that the electron temperature peak spreads out, ultimately 
creating the CF upstream. 

In the initial phase the ions are heated without appreciably affecting the electrons. Their 
dynamics can be approximated by that of a single fluid with half the pressure and whose 
sound speed is thus c^i = c^i/v^; we call this the ion- fluid. It is a flow with an effective 
Mach number Mj = \/2Mi and thermal conductivity from the ions alone: Pr = Pr^ = 0.33 
rather than Pr = 0.012 for the single-fluid system. Figure |H] shows the state of a two-fluid 
system with Mj = 2.4 and an ion-fluid at comparable times. Owing to its lower sound speed 
the ion-fluid solution is plotted against a coordinate rescaled to i-^ = ia/V2. The two-fluid 
system was initialized with gradient scale A = bin and the ion-fluid with A = bin. 

The spatial distributions of temperature, density and pressure are generally peaked 
at the center, s = 0. Figure [9] plots these central values versus time for this particular 
case. Mi = 2.4. This along with Figs. [7] and [8] illustrate and explain the four basic phases 
which compose the shock development in the two-fluid system. Except in cases with large A, 
discussed below, solutions at all Mach numbers pass through the same phases but at different 
times. 

The flrst phase is ballistic collapse during which the initial velocity gradient is steepened 
and density and pressures begin to rise. In the present example this occurs for t < T\y = 2.1. 
In the second phase (2 < t < 10), called ion heating, the ions are heated as they would be 
in a pure ion-fluid. The ion-fluid temperature will approach the steady value corresponding 
to a larger Mach number. Mi = \/2Mi. This is greater than, T2, so the ion temperature 
overshoots its final Rankine-Hugoniot value. At very large Mach number, T2 = 2T2, so the 
overshoot can be considerable. The left column in Fig. [7] consists of times during ion heating. 

In the third phase, front development (10 ^ t < 100), the ions transfer heat to the 
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Fig. 7. — Temperature of the ions (solid) and electrons (dashed) at times during the devel- 
opment of a shock. Inflow is at Mj = 2.4 with initial gradient scale A = 5£ji. The first eight 
panels are from successive times from t = I.Oth (top left) to t = SOOxji (center right). The 
horizontal axis is different in each column. The bottom right panel is the stationary solution 
from Fig. H] (Mi = 3.5). Grey curves are the time- dependent solution from the panel above. 
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Fig. 8. — Snapshots from the evolution of the two-fluid system (solid), the ion-fluid system 
(dashed) and single-fluid system (broken) at t = 4 (left) and t = 60 (right). The two- fluid 
and single-fluid systems have Mj = 2.4 while the ion-fluid has Mi — V^Mi = 3.40. The 
bottom panels show the temperature against a logarithmic length variable. The upper and 
lower solid curves are the ion and electron temperatures respectively. The center panels show 
the density at the same times. The top panels show total pressure (sohd), the single- fluid 
pressure (broken) or the ion pressure of the ion-fluid. The top axis shows Mm for the case 
Uei = 10^°cm-3 and Ti = 10^ K. 
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Fig. 9. — Time history of the central values for the Mj = 2.4 flow. The left panels show 
the temperature (bottom), density (middle) and pressure (top) vs. time; the bottom axis 
is in rescaled units and top axis in seconds for riei = 10^^ cm~^ and Ti = 10^ K. Solid 
curves are from the two-fluid solution; both electron temperature (lower curve) and ion 
temperature (upper curve) are plotted in the bottom row. Dashed and dash-dotted curves 
are the corresponding ion-fluid and single-fluid solutions. A grey line shows the peak ion 
temperature where that is different from the central value. The total pressure {pe + Pi) is 
plotted in the top row. A dotted line in the middle and top rows show isothermal, ballistic 
evolution which is singular at t = rt,. Symbols show the values at times where profiles are 
plotted in Fig. |8] (diamonds and squares) or in the right panels (crosses). The right panel 
show snapshots at t = 10 in the same format as Fig. [8] as well as the peak values (crosses) 
which are represented to the left. 
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electrons who transport it via thermal conduction to a growing CF. At the end of this phase 
the density and total pressure, Pe + Pi have achieved their final values, but the electrons 
are still colder than the ions at the center. The middle column in Fig. [7] consists of times 
during front development. In the final phase (t > 100) the ion cooling region develops in 
the center. The central curvature in the ion temperature reverses and the new pair of peaks 
move outward with the sub-shock. The central concavity descends to form a growing ion 
cooling region (see the right column in Fig. [7]). The decreasing central ion temperature 
meets the increasing central electron temperature at a level slightly below their equilibrium 
down-stream value, T2; they then increase together toward this final value {t > 500). 

The single-fiuid system evolves in a notably different manner than the two fluid system. 
During the ballistic collapse phase the density of all systems rise passively as 

. ^ 1 
n(0,t) 



1 - t/Tb 



shown as a dotted line in Fig. [91 This ends when pressure builds to a level sufficient to 
end the collapse. The collapse progresses farther in the single-fluid system because its large 
thermal conductivity makes it approximately isothermal, p ~ (1 — t/T\,)~^. In contrast to 
this, collapse in the pure-ion fluid, and thus in the two-fluid system, is better approximated 
as adiabatic p ~ (l-t/rb)~^/3. It is for this reason that the central density in the single-fluid 
system overshoots the Rankine-Hugoniot value far more than in the other two systems. 

Following the initial density overshoot, the central region of the single-fluid solution 
is forced to expand {t > 4). The expansion has a tendency to drive down the central 
temperature, and thus requires additional conduction from the shock to offset it. This 
competition creates a concave T proflle in the center, unrelated to ion cooling (obviously), 
which slows down the approach to flnal temperature. In the Mi = 2.4 case shown in Fig. [9] 
the flnal value is not achieved until after t = lOOOrji. The ion-collision time in the post-shock 
fluid is shorter that by a factor, (T2/Ti)^^/^ = 0.017, so the approach is truly slow. 



5.2. The Ballistic collapse 

The scenario above is altered in cases of extended ballistic collapse. The collapse time, 
Tb = A/Ml, is controlled by the initial gradient scale A. For very small values, A < 1, viscosity 
smoothes the initial gradient very rapidly and all subsequent evolution is unchanged. Larger 
values, on the other hand lead to a modifled scenario illustrated by the upper curves of Fig. 

For values of A large enough, the collapse phase can overwhelm the ion heating phase 
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Fig. 10. — Profiles of the ion (solid) and electron (dashed) temperatures at characteristic 
times for the flow with Mj = 2.4. Profiles are shown for fiows initialized with different 
gradient scales, arranged in rows with A increasing upward from A = 1 (bottom row) to 
A = 200 (top row). The vertical axis on each plot has tick marks at Ti and T2. 
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with consequences visible in the upper two rows of Fig. [TOl The slow collapse drives up both 
ion and electron temperatures adiabatically - they increase together. Electron conductivity 
then begins to generate a CF even before the ions have been heated. A shock ultimately 
forms away from the center, leaving the central ions unheated. These then cool as they 
subsequently expand, driving down the central ion temperature. This central depression 
superficially resembles the ion cooling region, but due to its different formation process it is 
a deeper concavity and the ion temperature can actually fall below the electron temperature 
(the upper right profile, A = 200, t = 300). Nevertheless, the shock conforms to the steady- 
state morphology and approaches the same steady profile. Its slower formation means the 
shock is slightly behind the location it would have had with smaller A. 



5.3. Motion of the shocks and fronts 

Shock tube solutions at different Mach numbers show the same development described 
above, but at different times. Figure E] shows that the extents of regions in the steady-state 
solutions scale strongly with Mi, and therefore with Mj. Figure [TT] shows how these regions 
move in time-dependent solutions. Each expands from zero to its final size. The higher Mach 
number cases have farther to go and thus develop more slowly. The ion cooling region is the 
largest and its development begins (symbols) substantially later for larger Mach numbers. 
The CF has attained its final size once its leading edge moves with the sub-shock — i.e. two 
curves in Fig. [11] have become parallel, as they are for t > SOOxji in the Mj = 3.6 solution 
(dashed curves). 

The different front development in the single-fluid system affects the motion of its CF 
as shown in the top panel of Fig. [TTJ Nearly instantaneous diffusion steals some of the 
post-shock pressure from the single-fluid sub-shock, causing it to lag behind its two-fluid 
counterpart. While the final shock speeds are identical for the two cases (from the Rankine- 
Hugoniot conditions) its higher pressure early on gives the two-fluid case a lead it never 
loses. The same instantaneous diffusion in the single-fluid system propels its CF ahead of 
its two-fluid counterpart at early times. This development appears more clearly diffusive, 
scF ~ y/t, than the two-fluid does. The steady-state CF width is, however, greater in the 
two-fluid case (see e.g. Fig. H]), so it eventually regains the lead. For the Mj = 5 case in the 
figure, it does this at t ~ SOOtji. Finally, it is clear that front development lasts far longer 
in the two- fluid case: the single- fluid curves are parallel at a much earlier point [t ~ SOOtji) 
indicating the CF has attained its final width. 
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Fig. 11. — Positions of the conduction front (CF) and sub-shock over time for shock-tube 
solutions. The bottom panel compares two-fluid solutions with different Mach numbers: 
Mi = 2.4 (solid), 3.6 (dashed) and 5.0 (broken). For each solution the upper curve is the CF 
and the lower is the sub-shock. Symbols on each pair show the beginning of the ion cooling 
phase. The top panel compares Mj = 5.0 solutions of the two- fluid problem (broken) with 
the single-fluid solution (solid). 



5.4. Validity of the fluid approximation 



During the evolution toward steady state, the flow develops gradients whose scale may 
challenge the assumptions underlying fluid approximation. For example, gradient length 
scales shorter than the mean free path violate the Chapman- Enskog derivation of classical, 
local heat flux. We demonstrated above that, at least in steady state, the electron temper- 
ature gradients are gradual enough to justify this approximation. We find the same to be 
true during th e time-evolution. This finding contradicts the hypotheses of rnany p revious 
investigations (IBrown et al.l 1 19791 : ISmith fc LilliequistI Il979l : iKarpen &: Devord 119871 ) which 
proposed modified treatments of the electron heat flux to handle the strong heating expected 
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in flares. Our findings are different because, while those investigators assumed electrons were 
heated directly, our electrons are heated through collisions with a much hotter ion popula- 
tion, which had themselves been heated in the shock. Since that coUisional coupling is far 
weaker than the conductivity, strong gradients never develop in the electron temperature. 

Since it is the ions that are first heated it is their temperature gradient that is most 
likely to violate the fluid approximation. Indeed, we find that Fi/Pj^^^^ becomes comparable 
to, and sometimes slightly greater than, unity during early phases of shock development. 
This is symptomatic of a more serious issue involving the viscosity itself. 

Viscosity is the ultimate source of heat but also contributes to decelerating the flow 
through the viscous stress. The pressure combines with the parallel co mponent of the viscous 



stress tensor, ass, to produce a diagonal element of the pressure tensor (ILifshitz fc Pitaevskii 



19811) 



Pss = Pi- CTss = Pi - ^fii-^ . (39) 

This is a second moment of the ion distribution function which can never be negative. This 
means that the velocity gradient must satisfy the condition 



4 /ij dv dv 
^ Pi ds ds 



Xi^ < I , (40) 



where the second expression involves rescaled variables. 

The velocity profile in a shock is principally compressive, dv/ds < 0, so condition ( HOl) 
is always satisfied. Even though we need not consider perpendicular momentum in our 
one- dimensional shock tube problem, there are diagonal perpendicular terms in the pressure 
tensor which must also be non-negative. The terms in the viscous stress, — cTss/2, will be 
positive in the case of compression and partially cancel out the pressure. Demanding non- 
negativity of this term completes the condition on the velocity gradient 

-2 < x^l^ < 1 , (41) 



where the upper bound is Eq. (140)) . 

Violation of criterion (jHJ) means one of the velocity second moments is negative. Since 
the squared velocity can never be negative this indicates that the distribution function itself 
is negative over a significant region in velocity space. While this is clearly unphysical, it is 
a common resul t of applying the Chap man- Enskog procedure beyond its range of validity 



(jCampbelllll984i ): the equilibrium Maxwellian has been supplemented by a negative pertur- 



bation so large the sum itself becomes negative. Nor does a positive value of the second 
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moment prove the distribution function was everywhere positive. Only for <^ 1 

can we be sure this is true, and that the classical, local transport formulae apply. 

Figure W2\ shows Xi9v/ds for all times and positions from a solution with Mj = 5 and 
A = 500£ji. It is small for t < rt = A/Mj = 100, and soon therafter violates criterion 
( BTj) . The violation is short-lived and confined to s < 300, but the region in the vicinity of 
the sub-shock continues to harbor values above half the limit; classical viscosity is probably 
inaccurate there. Outside the sub-shock itself, however, the condition is well satisfied and 
the fluid solution is probably valid. 
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Fig. 12. — A plot of the rescaled viscous stress, Xidv/ds, during the time development of a 
shock from Mj = 5 and A = bOOln. The quantity is plotted as a logarithmic grey-scale in the 
bottom panel against time (vertical axis) and position (horizontal axis). Solid contours are 
for values of —1 (outer) and —2 (inner). Dashed curves are the locations of the sub-shock 
(upper) and CF (lower). Horizontal dotted lines indicate three different times at which the 
viscous stress is plotted in the top panel (solid curves). Symbols show where the value crosses 
— 1 in each slice. The velocity itself is shown as a dashed curve at all three times, plotted 
against the right axis. Note that both axes are negative going upward, since both viscous 
stress and velocity are negative. 
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Shock solutions using methods beyond the fluid closures find very similar structure 
with slightly thicker sub-shock (ion shock) embedded inside the CF and ion cooling region. 
Steady sta te solutions have be e n found •with th e Fokker-P lanck equation in the steady-state 
method of iMott-Smithl (Il95ll ) (lGreywalllll975l : lAbelll975l ). Time-dependen t solutions have 



been found using fully kinetic treatment of ions along with fiuid electrons (jCasanova et al. 



199ll : IVidal et al.lll993l ). All such treatments find solutions very similar at modest Mach 
numbers (Mj ~ 2). Solutions at larger Mach numbers (Mj > 5) have sub-shocks 2-4 times 
thicker than fiuid solutions, but are otherwise extremely similar. We thus consider our 
fiuid solutions to be reasonable representations of the actual behavior in post-reconnection 
outfiows. 



6. The Slow shocks and conduction fronts in a reconnection jet 

The time-dependent shock solutions can be placed back into the context of Petschek 
reconnection to exhibit the structure of an outfiow jet . While the analysis has been time- 
dependent, we place the time-evolving solutions into a traditional steady-state model viewed 
from the end-on perspective of Fig. [TJ 



6.1. The jet geometry 

The configuration is structured by the rotational discontinuity (RD) which propagates 
along the inner field line, against the infiow, at the Alfven speed. Its location in the rescaled 
shock-tube variables, s and t, is 

SRD = (—-m) t = (.[^-m) t = ./^cos(A^/2)t , (42) 
\Csi J V V Po / V Po 

where is the angle between field lines across the initial current sheet (see Fig. [1]) and the 
final expression uses Eq. ([2]) for Mj and takes f2o,i = A6'/2. Until now we have assumed /3o 
to be small enough that sru was effectively infinite, or at least was ahead of the CF. 

Inside the RDs (s < sb.d) field lines have been advected downward at the vertical Alfven 
velocity, WAOy = f ao sin(A6'/2). They have been advected for At since reconnection at the 
X-point, giving them a vertical position 

y = VAOyAt = VAosm{A9/2)—t = J ^ sm{Ae/2) iat , (43) 

Csl V Po 

where t is the re-scaled time since reconnection (t = 0) and y increases downward. For 
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simplicity, we assume that the X-point is stationary, but it would be easy to add to Eq. 
an arbitrary X-point motion. 

The relation between the shock-tube coordinate s and the horizontal x coordinate nor- 
mal to the current sheet depends on aspects of reconnection we have not had to consider 
until now. The magnetic field n ormal to the cu rrent sheet, B^o, depends on the external 
structure in steady state models f|Petschek Il964h or the fast magnetosonic rarefaction wave 
in transient models (ILin &: Ledll994l : iNitta et al.ll2002l ). It depends on the reconnection rate 



(heretofore unspecified) and thereby determines the angle a = atan^B^Q/ Byo) the field lines 
make with the vertical in the end-on view. External solutions typically find this to be of 
order 10° or less. The value has no physical significance for our calculation so we adopt 
a = 10° for clarity of illustration. 

The field lines inside the RDs are horizontal and their length is parameterized by our 
rescaled shock-tube coordinate s. These horizontal field lines, with strength Bi, make an 
angle asin(Sa;i/i?i), with the line-of-sight (the z axis). The horizontal coordinate in the 
end-on view is therefore 



X 



Bxi 



Bxi B, 



il S 



ByO Bi 



where the final expression uses the fact that Bi = Bq and B^i — B^o 
[2l)-f H^ gives the half-angle the RDs make in the end-on view 



iis = tana sin(A6'/2)^iiS , (44) 

Combining expressions 



tanacos(A6'/2) 



(45) 



RD 



which is naturally smaller than a since the external field lines must intersect the RD. 



To cast the reconnection scenario in context we use the following values for the pre- 
reconnection plasma: Ueo = 2 x 10^°cm~^, Tq = 2 x 10^ K and Bq = 10 G. These choices 
dictate the parameters used in the rescaling, /3o = 0.028, in = 11.7 km and Th = 0.049 
sec. We set the external field angle at a = 10° and choose the current sheet angle, A^, 
corresponding to a particular shock-tube Mach number; for example, to obtain Mj = 5.0 
we use AO = 114°. Mapping the time-dependent Mi = 5 solution, with A = bOOin = 5.9 
Mm, onto the reconnection outfiow geometry yields the emission measure and individual 
temperature shown along the top row of Fig. [T3] f|T^ -d). 

The outfiow jet s differ in predictable ways from solutions of the full MHD equations 
previously reported (IForbes fc Priestlll983l : lYokoyama fc Shibatal 119971 ). The emission mea- 
sure (EM) is enhanced only within the center-most region confined by the sub-shocks. There 
are temperature enhancements outside this, representing the CPs. The profile of extends 
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Fig. 13. — Outflow jets synthesized from sliock-tube solutions at different conditions. Top 
row: = 114° (Mj = 5) sliowing (a) the regions, (b) the emission measure, (c) and (d) 
Tj. Bottom row = 95° (Mj = 3.6) showing (e) emission measure and (f) a combination 
of temperatures. Panels (g-h) are for = 77° (Mj = 2.4), from (g) two- fluid and (h) 
single-fluid solutions. 
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the farthest, but T, reaches the greatest value. The overshoot of Tj creates high-temperature 
ridges following the sub-shocks, between which is a cooler channel from the ion-cooling re- 
gion. Both EM and Tj are most enhanced within a short narrow region close to the X-point. 
This is a result of ballistic collapse due to the very large initial gradient- scale, A = 500£ji. 
All of these regions are, however, inside the RDs and therefore inside the outflow jet itself 
(dashed lines). This is the most notable difference from the simulations which were done for 
strictly anti-parallel reconnection, AO = 180°. 

The bottom row in Fig. [12] shows how the jets change at different Mach numbers and 
gradient scales. All use gradient A = lOin = 0.1 Mm, much smaller than the top row. 
As a result the sub-shocks connect to the X-point and the EM lacks an upper enhanced 
region. Instead the EM increases monotonically away from the X-point. The Tj ridge is 
progressively muted and the IC becomes progressively narrower moving from the largest 
reconnection angle {A9 = 114°, Fig. [T3li) to the smallest (A^^ = 77°, Fig. [T3^). At angles 
below this the region is so narrow that no structure is visible unless the horizontal axis is 
exaggerated — hence we do not show them. 

The two panels at the lower right (IT^-h) show the difference between electron temper- 
ature in a two-fiuid f|T3^) and single-fluid solution (11 3h). The main difference is that the 
slower rise in the former displaces the high-temperature region from X-point. Since imaging 
observations generally use radiative signatures sensitive to Te (line or continuum emission, 
for example), they would find a high-temperature "jet" separated from the X-point by dis- 
tances of 8-15 Mm (for Mj = 2.4 and 3.6 respectively). Single-fluid models fail to predict 
this since the electron temperature follows the ion temperature which rises much earlier. 



6.2. Extent of the conduction front 



A significant new feature of our solutions is that the conduction fronts all lie inside the 
RDs which d efine the outflow jet. This is notably different from simulations of anti-parallel 



reconnection, lYokoyama fc Shibatal (119971 ) in particular, where CFs extend outside the jets 
where they can drive evaporation. Our shock-tube treatment assumed, to the contrary, that 
all effects from the slow shock, including the CF, remained inside the RD. Our placement 
of the inflow boundary (the RD) at s — )■ +oo means that the shock-tube solution cannot 
itself violate the condition. It is, however, possible for the assumption to be violated once 
we commit to a particular choice of /3o in order to place the solution in a reconnection jet. 



Expression (j^2l) . giving the shock-tube coordinate of the RD, can be solved for (3o to 
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/3o = 1.2 ('-^ + MA . (46) 



yield 

B. = 1.2 ("-^ 

For a shock-tube solution at a given Mj we can find the position, scrit), of the CF over 
time, as shown in Fig. [TTl Setting srd = scf gives the condition for marginal validity of 
our shock-tube solutions. Using this in Eq. (H6ll yields the maximum value of /3o for which 
a particular solution satisfies the condition 

max(/3o) = 1.2 f^^^ + M,^ . (47) 



t 

Examples of this function are plotted on Fig. [TUfor solutions at different Mj and A = 10£ji 




Fig. 14. — The reconnection angle, max(A^) (top panel) or initial beta (bottom panel) at 
which the CF reaches the RD for shock-tube solutions at different Mach numbers. Expres- 
sions (H7j) and (HHj) are plotted against time (bottom axis) or against vertical coordinate in 
a steady-jet (top axis). Triangle marks a case with Mj = 5.0 discussed in the text. 

The shock-tube solutions are performed at specified inflow Mach number Mj. In this 
case the value of /3o corresponds to a reconnection angle through Eq. Using this in Eq. 
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f H7|) gives a maximum angle 



max(A0) = 4 sin ^ 



M2 ^^/^ 
max(/3o) 



(48) 



to which the internal CF assumption applies. The top panel of Fig. [T3] show the maximum 
angle for the different solutions. This shows that for Pq < 10^^ any reconnection angle 
< 120° will have CFs inside the RD, over most of the outflow jet. 

Since a CF decelerates over time, each of the curves in Fig. [13] trends upward. Time 
in the shock-tube solution is equivalent to vertical position in the outflow jet, according to 
Eq. f H3|) . It is thus possible for the CF to reach outside the RD in the upper portion of 
the jet (near the X-point) and return to the inside for the remainder of the jet. The section 
where scf > srd would involve a pre- heating of the material ahead of the RD whose effects 
we have not accounted for. For brief excursion we expect these effects to be minor since 
the propagation of the RD depends on density, which is far less affected by the CF than is 
temperature. Thus we can use the curves in Fig. [T3] to determine, at least approximately, 
the distance over which the CF extends outside the jet. 

For example, reconnection at /3o = 0.01 and an angle = 114° (dotted lines) produces 
an inflow of M, = 5. The CF position is outside the RD for t < 120 (triangle), corresponding 
to vertical coordinates y < 400\/5£ji- For Tq = 10^ and Ue = lO^'^cm"^, the CF extends 
outside the outflow over the first 5.2 Mm of the jet; beyond that position the CF lies inside 
the outflow jet. This does not occur in Fig. [T3k because that uses a solution with a large 
ballistic collapse, A = SOO^^i, for which the CF tends to lag. 

Other scenarios all follow the same general pattern. For a given reconnection angle, A^, 
lower values of f3o (i.e. larger Mi) will have larger CFs which extend outside the RD over a 
larger portion of the jet. Alternatively, for a given value of Pq larger reconnection angles will 
have more extensive CFs outside the RDs for more of the jet. The largest angle, A^ = 180°, 
produces a CF which i s everywhere outside the jet, as has been previo usly reported for anti- 



parallel reconnection ( iForbes et al.l Il989l : lYokoyama fc Shibatal 119971 ). (This limiting case 



cannot actually be treated using our shock-tube solutions. 



7. Discussion 

We have studied the structure of the slow magnetosonic shock in the outflow from 
Petschek reconnection. We have done this using two-fluid hydrodynamic equations, without 
including magnetic fields. This approximate treatment is possible for cases of reconnection 
between sufficiently skewed fields, A^^ < 150°. Such cases have been subject to less study 
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than the special case of anti-parallel field, = 180°. In skewed cases the primary role of 
the magnetic field is plasma acceleration at rotational discontinuities (RDs). This converts 
magnetic energy to bulk kinetic energy but without any heating. The heating occurs at a 
slow shock downstream of the RD; its primary source is the bulk kinetic energy from the 
accelerated fiow. It is a nearly parallel shock which we have treated as a purely hydrodynamic 
shock. 

Parallel shocks in ionized plasmas are best studied using two-fiuid equations. The den- 
sities and velocities of the electron and ion fluids are strongly coupled through electrostatic 
interactions, but their temperatures are coupled weakly though very slow collisions. This 
results in a structure more complex than what would be found from a single-fluid treatment, 
such as MHD, wherein electron and ion temperatures are assumed equal. The most signifl- 
cant feature of the two-fluid treatment is an extensive ion cooling region downstream of the 
ion deceleration layer (sub-shock). The slow development of the region, on the ion-electron 
collision time scale, will be reflected in the region of hot electrons within the outflow. Since 
most flare signatures, such as bremsstrahlung radiation, reflect the electron temperature, 
observations are better interpreted using two-fluid solutions than single-fluid (MHD) results. 
We have exhibited several examples of outflow jets synthesized from two-fluid hydrodynamic 
solutions. 

The most signiflcant difference between our solution and previous studies comes from the 
skewed flelds rather than from the two-fluid treatment. The separation between RD and slow 
shock, which occurs in all skewed flelds, means the conduction front is often located inside the 
outflow jet. In this case the chromospheric evapora tion will not precede the reconnection, as 



it does in the anti-parallel cas es previously studied (IForbes et al. 



19891 : lYokoyama fc Shibata 



19971 : ISeaton fc Forbes! l2009l ). Evaporation has, however, been postulated to enh ance the 



emission measure of high temperature material to levels actually observed in flares (ITsuneta 



19961 ). We conclude that this assertion should be made in conjunction with an estimate of 
the strength of the current-aligned fleld, sometimes called the guide fleld, at the reconnection 
site. 

By using purely hydrodynamic equations we have entirely neglected the energetic contri- 
bution from the magnetic fleld decrease across the slow magnetosonic shock. With constant 
fleld strength the source o f kinetic energy in our model must be decreasing fleld line length 



see 



Longcope et al.l l2009l . for an analysis of this). An initially bent field line relaxes to a 
straight line which is shorter by a factor cos(A6'/2), ignoring the current sheet thickness. 
The magnetic energy decreases by this same fraction, provided the field strength does not 
change. Had we accounted for the field decrease at the SS the fraction could be smaller still, 
so ours is a conservative estimate of the energy release. On the other hand, the antiparallel 
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limit {A9 = 180°) predicts a 100% magnetic energy release from shortening alone, to which 
the field strength decrease can add nothing. 

Figure [2] shows that for angles A^^ > 120° the post-shock field strength will decrease by 
more than one-third; this decrease is missing from our calculations. The top pan el shows 



that properly accounting for this decrease (solid curve, for the full MHD solution of ISoward 



19821 ) yields a post-shock temperature lower than that in a straight shock tube (dashed) 
which ignores it (albeit only slightly lower). The reason that greater energy release results 
in lower temperature is that the fiux tube is wider where it is weaker. The extra magnetic 
energy released by weakening is more than fully consumed by work done expanding the tube. 
The very similar post-shock temperatures in the two cases, in spite of the expansion in the 
first, leads us to believe that our shock-tube hydrodynamic model is adequate for studying 
thermal effects. Neglecting magnetic energy release appears, paradoxically, to produce a 
slight overestimate of the temperature in the CF. 

Our two-fluid modeling uses viscosity and thermal conductivity based on Coulomb col- 
lisions. Notably absent is coUisional resistivity which has played a central role in previous 
reconnection models. Resistivity would have played its most most significant role generating 
the reconnection electric field within the diffusion region. Our analysis has assumed this 
event to have occurred before t = 0. We have neglected direct energetic contributions from 
the resistive dissipation accompanying the electric field (i.e. Ohmic heating). We neglect the 
contribution from this very small region in order to focus on the heating at the slow shocks. 

Nor would resistivity have significantly affected the shocks, had it been included. Since 
the magnetic field bends only slightly at the slow shocks (see Fig. [2]), they carry little current 
and therefore account for little resistive dissipation. In our shock tube treatment there is 
none at all. More significant current will occur at the RDs, where the magnetic fields 
are bent following reconnection. Including resistivity in the RD dynamics would diffuse 
these structures, replacing sharp bends with smoother corners. In this case they are time- 
dependent structures rather than conventional steady shocks. Diffusion would broaden the 
RDs, diminishing the direct Ohmic heating within, but would leave their global structure 
unaffected. Thus we reach the conclusion that resistive dissipation makes a minor energetic 
contribution in skewed reconnection. 

We find the use of coUisional transport to be warranted for electron thermal conduction. 
The rate of electron heating is low compared to the conductivity, so electron temperature 
gradients are always shallow. The ion heating, directly from viscous dissipation, is far more 
rapid and leads to gradients at the limit of coUisional ( Chapman- Enskog) treatment. This is 
related to the even more significant issue of the steep velocity gradients from the ion pressure 
and viscous stress. These are also at the margin of fiuid treatment; we found instances where 
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the Chapman-Enskog treatment was unphysical, although only temporarily. Thus we expect 
non-fluid effects to be significant for ion dynamics within the sub-shock. Previous studies 
have determined that such a treatment leads to a sub-shock thicker than that from fluid 



closures, by a factor in the range of two to four (lGreywalllll975t IVidal et al.lll993l ). Aside 
from this thickened inner region, we expect our solutions to be approximately valid for cases 
of skew angle M < 150°. 

Some recent modeling of magnetic reconnection has f ocused on electric fie l d from col 
lisionless effects such as electron inertia or the Hall effect (IBiskamp et al.l llQQTt IShay et al. 
19991 : iBirn et al.ll200ll ) . While mo st are two-dimensional stud i es of anti-parallel reconnection, 
some have included a guide field ( iPritchett fc Coronitil l2004l : iHesse et al.l |2004J : iDrake et al. 
20061 ). The coUisionless mechanisms generate an electric fields localized to a region of col- 
lisionless scale, such as di = c/upi. The coUisional scale of our solution is larger than this 
by 

in Vth 



di c 

where A = n\\) ^ 1 is the plasma parameter and Ad is the Debeye length (lNicholsonlll983l ). 



A 
h^A 



(49) 



For the parameters we have used, rigi = 10 cm~ and Ti = 10 K, the ion skin depth is 
di = Q X 10~^£ji- This means the mean-free path is effectively infinite compared to dj, and 
the plasma appears coUisionless on this scale. Conversely, the coUisionless effects occur on 
scales far smaller than all the structures in the outflow jets, including the slow magnetosonic 
shocks. 

CoUisionless effects would become significant were flux transfer to occur within a region 
of coUisionless scale. This would enter our post-reconnection analysis as an initial gradient 
scale A ^ (ij ^ We have found that all initial conditions with A < tn converge to those 

— 1/2 

we have analyzed within several ion collision times. This is within a distance Ay ~ ^ii/^o 
of the X-point itself, beyond which the solution would not differ from a purely collisional 
solution where A ~ fji. Thus we expect the large-scale appearance, including the slow shocks, 
to be similar even when coUisionless processes are responsible for the flux transfer. 

The initial relaxation of the velocity gradient in our fluid model is through viscous 
diffusion. This would not be a valid approach were the initial gradient scale very small: 
A ^ in. Ion-ion collisions would be ineffective over such small scales and the counter- 
directed parallel flows from the two RDs would interpenetrate, forming a central region of 
counter-streaming ion beams. Such counter-str eaming ions have b een found in simulations 
(jPrake et al.ll2009h and spacecraft observations (iGosling et al.ll2005l ) when reconnection flux- 
transfer occurs on sub-coUisional scales. Calculations suggest that plasma instabilit ies will 
devel op to thermalize the ion distribution and form a hot, stationary, central region ( iParker 
196ll ). Since any such mechanism must satisfy the same conservation laws, the final structure 
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will be identical to that found using two-fluid equations. Only the thickness of the sub-shock 
will differ owing to the modiflcation to the ion momentum diffusion (i.e. viscosity). 

We believe our analysis is fairly robust and broadly applicable. It models the dynamics 
following a sudden change to a fleld line's topology which was assumed to occur within 
a small region - this is the flux transfer. This topological change creates a mechanical 
disequilibrium of the fleld line causing it to shorten on the Alfvenic time scale. The rapid 
shortening will drive compressive flows far faster than the sound speed (provided (3 <^ 1). 
It is the collision of these super-sonic compressive flows which generate the slow shocks and 
associated conduction fronts. While we have studied the dynamics in a steady state, two- 
dimensional model, the same scenario has been identifled in t ransient, three-dimensional 



reconnection (ILongcope et al.l l2009l : iGuidoni fc Longcopd |2010| ) . Even in the absence of a 



steady, coherent outflow jet, we expect Alfvenic retraction, slow shocks and conduction fronts 
to be inevitable features of fast magnetic reconnection. 
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